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THE  ART  OF 

BOND  GRAPH  CONSTRUCTION  FOR  TRANSDUCERS 


Bulletin  #1  -  r.ongitudinal  Vibrator  (A  2-Port  Model) 

(Ref.:  F.  Hunt  “Electroacoustics”  and  S.  Hanish  “Treatise  on  Acoustic  Radiation”  Vol.  II NRL) 

The  longitudinal  vibrator  can  conveniently  serve  as  a  prototype  in  the  art  of  bond  graph  construction. 

Suppose  first  we  have  a  moderately  long  thin  elastic  bar  of  constant  cross  sectional  area  and  assume  it  is 
in  steady  state  vibration  at  a  fixed  frequency.  We  propose  to  construct  a  bond  graph  model  of  this  bar.  Our 
procedure  is  outlined  in  a  sequence  of  steps. 

(1)  As  an  initial  step  we  perform  a  dynamic  analysis  of  the  motion  of  the  bar.  But  immediately  we  note  a 
different  approach:  in  conventional  analysis  we  formulate  an  equation  of  velocity  resulting  from  applied  force; 
in  a  different  approach,  because  a  bond  graph  pictures  a  flow  of  power,  we  seek  to  find  the  velocity  at  every 
point  in  the  bar,  as  well  as  the  force  at  every  point,  in  terms  of  the  force  and  velocity  at  one  end  of  the  bar.  Let 
us  set  aside  the  physical  causes  of  this  force  and  velocity  at  the  input  point  (that  is,  bar  end).  Because  the  bar 
is  an  elastic  continuum  the  force  and  velocity  at  the  second  end  are  easily  obtained.  Analysis  therefore  leads  to 
two  equations  relating  forces  Fj,  F^  and  velocities  Vj,  V^. 

F=F,(F,,V.) 

V,=V,(F,,V,) 

In  most  common  applications’  we  find  the  functional  relations  to  be  linear,  and  the  power  variables  decoupled, 
i.e. 


F2=AF,  +  BV, 

V,  =  CF,  +  DV, 

Here  A,B,C,D  apply  only  to  the  terminii  of  the  bar.  By  so  doing  we  have  converted  the  bar  description  from 
that  of  a  continuum  to  a  system  of  2  degrees-of-freedom,  namely  the  velocities  of  the  ends  of  the  bar.  On 
dimensional  considerations  it  is  seen  that  A,B,C,D  are  impedances,  or  admittances,  or  ratios  of  impedance’s. 

The  above  set  of  equations  in  F^,  V^,  obtained  by  dynamic  analysis,  has  no  relation  to  any  preconceived 
model,  be  it  electrical  analogies  or  bond  graphs,  or  signal  graphs,  or  chose  what  one  will. 

(2)  The  second  step  in  the  procedure  introduces  modeling  by  bond  graphs.  To  begin,  one  identifies  each 
one  of  two  degrees  of-freedom  with  a  “one-junction”  (written  as,  -1-),  and  because  these  are  coupled,  the 
coupling  (by  bond  graph  rules)  is  a  “zero-junction”  (written  as  -0-).  The  primitive  graph  is  then. 
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The  digits  identify  the  power  bonds,  which  are  paths  of  power  flow.  Bonds  #1,  #7  are  “inputs,”  that  is  “ports.” 
Bonds  #2,  #4,  #6  are  reactive  impedances  or  admittances  associated  with  a  lumped  parameter  2  degree-of- 
freedom  system.  Bonds  #3,  #5  are  interior  coupling  bonds.  Completing  a  more  elaborate  graph,  one  has. 


Z2  Y4  Z6 

4  4  4 


F, 


V,  F3 
Fb 


V,  V. 


0 


v,  F* 

F3 


V, 


V2 


The  arrow  directions  are  arbitrary,  but  it  is  advantageous  to  assign  them  “all  right  pointing,”  “all  left  pointing,” 
or  “symmetrical,”  (as  above),  Z,,  are  mechanical  impendances,  and  Y4  is  a  mechanical  admittance. 

The  above  graph  now  permits  writing  of  “graph  equations”  in  accordance  with  the  bond  graph  rules 
(1)  that  1 -junctions  represent  the  operation  of  force  summation,  and  (2)  the  0-junctions  represent  the  operation 
of  velocity  summation.  The  “graph  equations”  are  then: 

(a)  F,  =  Z,V,  +  F3;  (b)  V.  =  Y,F,3  -  V^;  (c)  F,  =  Z^V^  +  F3 
Here  the  coupling  force  Fj^  is  the  same  as  as  F3. 

By  construction,  the  model  has  now  three  unknowns,  namely  Z,,  Y^,  Z^. 

(3)  The  third,  and  crucial  step  in  the  construction  of  the  bond  graph  model  is  to  rearrange  the  “graph 
equations”  so  that  they  appear  as  the  set, 

F,  =  GF,  +  HV,;  =  JF,  +  KV, 

in  which  G,  H,  J,  K  are  algebraic  functions  of  Z^,  Y^,  Zg,  and  are  therefore  unknown.  To  find  them  explicitly 
we  return  to  the  dynamic  equations  obtained  by  analysis,  and  equate  comparable  coefficients: 

(a)G  =  A;  (b)H  =  B;(c)J  =  C;(d)K  =  D 

Because  A,  B,  C,  D  are  known,  this  set  allows  Z^,  Y^,  Zg  to  be  known  function  of  one,  two,  three  or  all  four  of 
them. 

The  model,  as  a  first  approximation  (i.e.  2  degrees-of-freedom)  is  now  completed.  It  is  seen  to  be  a  “one- 
zero-one”  structure,  that  is,  a  conventional  “Tee”  in  electric  circuit  theory  which,  in  itself,  is  an  approximation 
to  a  short  transmission  line.  However  it  is  to  be  emphasized  that  the  bond  graph  model  is  independent  of 
electric  circuit,  or  electric  transmission  line  theory. 

(4)  The  fourth  step  in  the  procedure  is  to  introduce  modal  description  of  vibration  of  the  bar  by  expanding 
Z2,  Y4,  Z6  in  resonant  or  anti-resonant  modes.  This  is  possible  because  Z^,  Y^,  Zg  are  transcendental  functions 
of  frequency.  The  type  of  expansion  is  optional.  We  choose  here  to  expand  Z^,  Zg  in  resonant  modes,  and  Y^  in 
antiresonant  modes.  The  bond  graph  symbols  for  expansions  are  “fans  of  bonds.” 


Ri 


I2  C2R3 

1 V. 


'O' 


I^C3 


_Fj 

Vi  V, 


C2  G2 

G.  n  \/  r, 

\  /r  O' 

C,— 0 
F„ 


V, 


C3  G3 


0- 


c,  R, 
\/ 

"vN 


O' 


-1^^2 

Vb 


2 


Here  v„,  v^,  ,  are  modal  velocities;  F^,  Fp,  F^ . . .  are  modal  forces;  C^,  Cp,  . . .  are  modal  compli¬ 
ances;  R„,  . . .  are  modal  resistances.  Similarly,  C,  G,  r  are  inverses  of  mass,  resistance  and  compliance. 

All  modal  parameters  (I,  C,  R,  C,  G,  r)  are  known  from  the  results  of  step  #3  above. 

The  bond  graph  of  the  mechanical  motion  of  the  bar  is  now  completed. 


Explanation 

(1) 


(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 


Observations 

(a) 


the  lumped  parameter  model  is  a  cascade  of  masses  and  springs. 

t  h  ^ _ t'  b  ^ ^ — nnr\ — 

Vj  is  coupled  to  via  a  cascade  of  model  masses  and  modal  springs,  the  magnitudes  of 
which  are  frequency  dependent. 

at  any  frequency,  velocity  Vj  is  a  sum  of  modal  velocities  Vp(©),  V^(©),  etc.  which  are 
determined  by  modal  Z^((o),  Zp((0),  etc. 

at  modal  resonant  frequencies  to^,  (O^, . . .  etc.,  the  Z„((0^),  Zp(a)5),  etc.  are  minimized  so  that 
V  or  V- . . ,  etc.  is  maximized. 

a  p 

at  any  frequency  (O  the  force  FjCco)  which  couples  Vj  to  is  a  sum  of  modal  forces  F^((n), 
Fp(©),  F^(©),  etc. 

at  modal  antiresonant  frequencies  ©j,  etc.,  the ¥^^^(©'^1,  Yp(©2),  etc.,  are  minimized  so  that 
F„(©^),  Fp(©^),  etc.,  are  maximized. 

the  assignment  of  directions  of  power  flow  at  any  frequency  ©  is  purely  arbitrary,  but  consis¬ 
tent  choices  can  be  made. 

the  equations  of  motion  can  be  written  directly  by  summing  forces  at  1-junctions  and/or 
summary  velocities  at  0-junctions. 


If  the  bar  is  driven  by  an  external  source  of  power  the  model  then  contains  a  source  bond. 
Inspection  of  the  physical  requirements  shows  at  what  point  in  the  graph  it  is  to  be  attached. 
For  example,  if  the  bar  itself  is  piezoelectric  the  driving  force  in  a  2  degree-of-freedom  model 
shares  the  same  velocity  difference  as  the  admittance  Y^.  Its  representation  is  therefore. 


S— 1— Y, 
IV, -V, 
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(b)  Since  the  bar  is  an  elastic  continuum  more  points  of  velocity  can  be  used  in  the  dynamic 
analysis  as  needed.  For  example,  if  a  compliance  load  is  attached,  x  units  of  distance  from  the 
end  its  velocity  can  be  represented  (in  approximation)  as  the  third  degree-of-freedom  in  a 
3  degree-of-freedom  system,  i.e. 


V, 


(c)  Since  an  elastic  continuum  may  vibrate  in  several  different  modes — longitudinal,  shear,  twist, 
bulge,  radial  etc. — it  is  a  common  occurrence  of  a  vibratory  motion  that  two  or  more  of  these 
modes  exist  together  over  particular  frequency  ranges.  Cross  coupling  between  them  is  then 
to  be  expected.  For  example,  suppose  a  vibration  is  simultaneously  longitudinal,  shear,  and 
radial  and  assume  the  modeling  makes  the  following  assignments: 


longitudinal  mode:  2-d.o.f.  (Vi,,  Vy 
shear  mode:  2-d.o.t  (V^,,  V^j) 
radial  mode:  1-d.o.f.  (V„) 


(Note:  d.o.f.  means 
‘degree  of  freedom’) 


Assume  further  that  analysis  shows  cross  couplings  of  compliance,  inertance,  and  resistance 
between  respective  velocities.  An  illustrative  graph  may  have  the  following  appearance: 
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Cross  couplings  modes  are  therefore  easily  made  visible  and  tractable 
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An  Alternative  Approach 
A  3-Port  Model  of  a  Longitudinal  Vibrator 

(1)  Asa  first  step  we  anticipate  a  model  of  the  bar  to  have  the  character  of  a  multiport  representation.  A 
port  by  definition  is  a  location  on  the  surface  of  the  bar  where  power  may  enter  or  exit  during  motion  of  the 
bar.  In  bond  graph  theiry  a  port  is  pictured  as  a  short  straight  line  terminating  at  one  end  in  a  junction  of  the 
bond  graph  and  directed  at  the  other  end  away  from  the  graph. 

Inspection  of  the  bar  shows  the  existence  of  several  possible  ports.  We  chose  three;  one  at  each  end  of  the 
bar,  and  one  on  the  flat  surface  yj,.  Each  port  is  identified  by  its  power  variables  of  effort  and  flow  properly 
subscripted. 

(2)  A  second  step  is  to  perform  a  dynamic  analysis  of  the  motion  of  the  bar  for  the  condition  of  steady 
state.  The  objective  of  this  analysis  is  to  associate  with  each  port  an  explicit  mathematical  form  of  its  input 
impedance  or  admittance.  This  procedure  leads  to  a  multiport  solution  of  the  motion  of  the  bar. 

(3)  The  next  step  is  to  construct  an  elementary  bondgraph  of  the  vibrating  bar.  Each  port  (in  the  form  of  a 
short  straight  line)  is  attached  to  a  1-junction,  that  is  to  an  “effort  divider.”  Since  the  motion  at  the  ends  are 
coupled  to  each  other,  their  ports  (here  designated  as  A,  B)  are  connected  via  aO-junction,  that  is,  to  an  “effort 
divider.”  The  flow  of  the  third  port,  here  designated  as  port  C,  is  coupled  to  ports  A,  B  only  through  the 
already  established  0-junction.  The  primitive  bond  graph  can  now  readily  be  constructed; 


PortB^  ^  ^  Port  A 

Here,  the  symbol  is  the  input  self  admittance  of  port  B,  and  of  port  A.  The  symbol  h,  is  the  input  self 

admittance  of  port  C. 

(4)  Since  the  bar  in  question  is  mechanical,  its  power  variables  are  force  F  and  velocity  V.  E^r  various 
reasons  the  makers  of  models  have  found  it  advantageous  to  choose  F  to  be  “flow,”  and  V  to  be  “effort.”  A 
more  complete  bond  graph  incorporates  these  choices; 


PortC 


h, 

]v. 


5* 


Vb 


rw 


112 
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PortB 


0- 


Port  A 


(Note  that  arrow  directions  are  arbitrary,  but  conform  to  convention.) 
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(5)  We  next  formulate  the  constitutive  equations  relating  force  to  velocity  at  each  port.  For  the  choice  of 
effort  and  flow  variables  made  earlier  these  are 

(l)v  =  h  F;(2)V  =h,F-(3)v  =h,F 

Za^^'b  2b''w  Iw 

It  is  seen  that  h,,  h,  are  mechanical  admittances. 

(6)  Explicit  forms  of  hp  h,  are  available  in  the  literature  (Ref.  1).  Our  next  step  is  to  incorporate  them  in 
the  constitutive  equations 


a)v,= 


hpF, 

F  a 


J  tan 


(2)v,= 


Ml. 

j  tan^w/^j 


(3)  v^=j  hp  sin 


h  c2— i- 

^  Vd  pc  ’  P 

E=  modulus  of  electricity  |units: 
p  =  mass  of  the  bar  per  unit  volume  (units:  N  secVm'* 


(7)  The  bond  graph  model  is  now  completed.  For  the  choices  of  power  variables  it  is  seen  to  have  the 
form  of  a  “tee”  (i.e.  —  1  —  0  —  1  — ) 

(Ref.  W.R  Mason,  “Physical  Acoustics,”  Academic  Press,  New  York,  1964.) 

Example  (a  3-port  model  of  a  quarter-wave  resonator). 

If  velocity  is  made  zero  by  fixing  end  A  of  the  bar  in  a  wall,  the  graph  constructed  above  takes  on  the 
form. 


PortB 


p  Port  C 


Since  the  0-junction  now  has  no  exterior  bond  one  sets  F^  =  F^.  The  reduced  graph  is  now. 


h, 


Port  B  p  Port  C 

1  W 


The  admittances  h,,  h^  are  now  in  series,  each  having  the  flow  F^.  They  can  be  added  to  form  h, 
h=h.+h,=jh,si„(:^) .  _i  =jh,u„(:^) 
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Setting  V  =  (note  that  these  v’s  are  ‘across’  variables,  that  is  ‘efforts.’)  the  graph  is  simplified  to, 

h 

|v  v  =  hF^ 

Port  B  n  Port  C 
Pw 

Since  h  is  a  transcendental  function  it  can  be  expanded  in  modes.  We  choose  antiresonant  modes,  namely 
modes  wherein  at  any  as  a  sequence  of  antiresonant  frequencies  the  ‘across’  variable  (here  velocity  v)  is 
maximized  and  the  ‘through’  variable  (FJ  is  minimized.  The  graph  then  becomes: 


C(W,) 


=  a  quarter- wave  longitudinal 
resonator 


Other  examples  may  be  found  in  the  cited  references. 

Tauchspulenmikrophone 

Bulletin  #2  -  Electrodvnamic  Microphone  (Cascade  Model) 

In  the  book  “Grundlagen  Der  Sechnischen  Akustik”  by  W.  Reichardt,*  there  appears  on  page  305  a 
description  of  a  Tauchspulenmikrophone  (plunger  spool  microphone  =  electrodynamic  microphone).  Its 
essential  parts  are  reproduced  below. 


N,h; 


Sph 


Electrodynamic 
microphone 

Cross  section  of  electrodynamic 
microphone 

Here  n,  m,  h  are  lumped  compliance,  mass,  and  resistance  respectively. 


*  Akademische  Verlagsgesellshaf  Geest  Portig  K.  -  G.  Leipzig  1968 
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Acoustic  Description 

This  device  is  an  acoustic  pressure  receiver.  Sound  pressure  incident  from  the  left  passes  through  screen  b 
and  protective  grid  w  to  reach  the  flexible  diaphragm  a  to  which  a  spool  of  wire  is  attached.  Directly  behind 
the  diaphragm  are  a  series  of  cavities  and  orifices  which  perform  the  function  of  smoothing  the  frequency 
response.  The  spool  of  wire  is  transverse  to  a  static  (d.c.)  magnetic  field  created  by  permanent  magnet  c.  In 
the  rear  cavity  is  an  equalizer  tube  whose  orifice  f  admits  “atmospheric  pressure”  via  tube  feed  g. 

Bond  Graph  Model 

A.  Construction  of  a  bond  graph  model  of  this  device  requires  an  understanding  of  the  acoustic 
performance  of  cavities  and  orifices.  To  illustrate  this  performance  we  propose  to  make  a 
bondgraph  model  of  the  following  simple  cavity  and  orifice  structure. 


orifice 


C 


p2 

q2 


We  note: 

(1)  the  flow  q  is  hydrodynamic  compressible  flow,  q,  q^ 

(2)  the  orifice  length  1  is  “short”,  meaning  it  has  constant  flow  qj 

(3)  all  variables  pj,  p^,  q,,  q^  are  sinusoidal  steady  state. 

(4)  a  flow  resistance  R  is  assumed  to  be  present  in  the  orifice. 

(5)  all  pressures  are  relative  to  atmospheric  pressure. 

The  corresponding  bond  graph  is  easily  constructed: 


Explanation 

(a)  the  1-junction  (velocity  qj)  is  a  “pressure  divider”,  i.e.  the  input  pressure  p,  is  “divided”  by 
the  orifice:  part  p,  accelerates  an  air  mass  m  (units  NsVm’),  and  part  p^  (Ns7ra*)  overcomes 
the  orifice  resistance.  After  subtraction  of  Pj  and  p^  the  pressure  p^in  the  chamber  remains. 

(b)  the  0-junction  is  a  “volume  flow  divider”  i.e.  the  input  flow  qj  in  chamber  #2  is  partly  di¬ 
verted  in  amount  qj-q^to  be  associated  with  the  compression/expansion  of  the  chamber 
volume  acting  as  a  “hydrodynamic  spring”  of  compliance  C,  leaving  exit  flow  q^  as  remain¬ 
der.  The  compliance  C  equals  1/K  where  K  (^/m’)  is  the  “spring  constant.” 

(c)  the  1 -junction  (velocity  q^)  divides  p^  into  p^^  and  p^. 

(d)  the  directions  of  power  flow  through  the  junctions  reverses  in  the  second  half  of  the  sinusoi- 
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dal  cycle,  but  not  the  directions  of  power  flow  in  load  elements  1,  R,  C. 

(e)  the  assignment  of  causality  allows  the  formation  of  the  constitutive  equations  for  1,  R,  C: 

(1)  Qj  =  pjdt  (2)  pj^  =  Rq^^  (3)  Pj,  =  ^  J  q^dt 

B.  We  now  proceed  with  the  construction  of  the  bond  graph  model.  The  component  parts  are, 

(a)  (forward)  source  (d)  equalizing  tube 

(b)  diaphragm  (e)  (rear)  source 

(c)  (smoothing)  orifices  and  cavities 

(a)  (forward)  source 

Since  the  microphone  is  a  pressure  sensor  we  designate  the  exterior  driver  to  be  an  “effort 
source”  Se.  Because  of  the  presence  of  screen  and  grill  we  take  the  source  impedance  to  be 
the  acoustic  resistanceR^.  To  construct  the  bond  graph  of  this  source,  the  rule  is  invoked  that 
an  effort  source  has  the  operational  significance  of  a  0-junction,  and  therefore  must  be 
followed  by  a  1-junction.  The  model  is  therefore. 


Se 


R^  has  been  added  to  account  for  radiation  resistance,  if  not  negligible. 

(b)  diaphragm 

The  diaphragm  is  a  mechanical  component  characterized  by  the  power  variables  force  F; 
velocity  It  is  assumed  in  the  first  approximation  to  have  rigid  piston  motion,  featuring 
inertance  Id,  capacitance  C^,,  a  resistance  R^.  In  operation  the  pressure  difference  between  the 
opposing  sides  of  this  piston  excite  it  into  motion.  The  bond  graph  model  is  therefore  an 
“effort  divider,”  that  is ,  a  1-junction: 


(c )  (smoothing)  orifices  and  cavities.  For  identification  of  chambers  see  subscripts  on  cross- 
section.  The  first  cavity  behind  the  diaphragm  is  designated  by  the  symbol  (compliance 
#2).  Volume  velocity  q^  enters  and  velocity  qp  exits.  The  bond  graph  model  is  therefore  a 
“velocity  divider,”  i.e.  a  0-junction, 

Q 

Pckz 
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(d)  the  first  orifice  connects  chamber  #2  to  chamber  #3.  Its  volume  velocity  is  qp.  Pressure 
differences  between  the  two  sides  of  the  orifice  (tube)  accelerate  a  mass  m^  and  overcome  a 
resistance  R^.  The  bond  graph  model  is  therefore  an  “effort  divider”  i.e.  a  1 -junction: 


m3  R3 

i’s  iid 


%  qp 


(e)  the  second  chamber  #3  has  input  volume  velocity  Qp  and  exit  velocity  q^.  Its  model  is  there¬ 
fore  a  0-junction: 


Cj 


fl 

‘Ip 


‘1r-‘1v 


(f)  the  second  orifice  connects  chambers  #3  and  #4.  Its  graph  is, 

I  R 

J2 

(g)  the  third  chamber  #4  has  input  volume  velocity  q^  and  exit  velocity  q^  (=  equalizer-tube 
velocity),  the  graph  of  which  is, 

C4 

Psl 


(h)  the  third  “orifice”  is  the  entrance  and  length  the  equalizer-tube  in  which  the  velocity  by 
approximation  is  q^.  Its  graph  is. 


(i)  (rear)  source 

The  rear  source  is  the  external  acoustic  pressure  impinging  on  the  surface  aperture  of  the 
equalizer  tube.  Its  graph  is  structurally  the  same  as  the  forwarded  source  but  has  different 
magnitudes  of  the  power  variables  and  resistive  loads. 


Assembly  of  Components  (Cascade  Model) 

In  assembling  the  components  of  a  bond  graph  one  follows  the  simple  rule  that  bonds  having  the  same 
power  variables  coalesce  into  a  single  bond. 
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Assemblies  can  be  made  to  various  degrees  of  approximation.  We  sketch  first  an  assembly  which  neglects 
area  transformations  in  the  path  of  flow. 


|— forward  source  —\ - diaphragm - 1— chamber  #2— |— chamte*#S"#2  1 


C  I 


R 

/ 

Pr 


nif 

I 


Rf 

R  R 


Ps 


% 


-0- 


Ps 


(qv-q*) 


-0- 


p; 


<k 


q,'=qs 

rear 

source 


•s; 


In  a  second  approximation  area  transformation  can  be  modeled  by  use  of  simple  numerical  ratio  trans¬ 
formers. 

Rlectrical  Output 

Electrical  output  of  the  microphone  is  accomplished  by  Faraday  induction.  In  this  process  the  axial 
velocity  V„  of  the  diaphragm  is  imparted  to  the  “spool”  of  wire  so  that  the  wire  strands  move  normally  to  the 
(d.c.)  magWc  flux  B  which  is  generated  radially  by  a  permanent  magnet.  By  induction  an  e.m.f.  appears 
across  the  terminal  of  the  spool  of  wire,  length  1,  of  amount, 

Eo.c  =  (VxB)Jl 

For  wire  resistance  the  wire  electrical  current  is  ig=E^/R^.  Since  E^  is  operationally  an  “effort”  it  has 
the  character  of  a  negative  source  (i.e.  power  output  rather  than  power  input).  Its  bond  graph  is  then 


E.:S. 


E, 


E„ 


The  port  Fg,V^  is  bonded  to  the  1-junction  labeled  q^. 

Conclusion 

The  “plunger-spool”  microphone  has  been  modeled  above  as  a  lumped  parameter  device  with  all  compo¬ 
nents  in  cascade.  In  practice  the  diaphragm  will  vibrate  in  multimode  patterns  of  surface  displacement.  These 
can  be  modeled  as  a  “fan  of  resonant  modes”  described  in  Ref.  2. 
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Bulletin  #3  -  Transducers.  General  Principles 


In  bond  graph  theory  transducers  are  modeled  as  multiport  systems  in  which  the  power  variables  of  bonds 
are  interrelated  by  coefficients  arranged  in  matrices.  For  a  system  of  n  ports,  in  which  the  power  variables  of 
the  I’th  port  are  ei,  fi  (effort,  flow),  the  interrelations  are  expressed  by  these  matrix  equations: 

(a)  to  first  (or  linear)  order  in  ej,  fj: 


n  n 

Cj  =  X  Aij  Q  +  X  Bij  ej,  i  =  1,2, . . .  n 

j=l 

n  n 

fj  =  X  CiJ  Q  +  X  Dij  ej,  i  =  1,2, . . .  n 

j=i 

(b)  to  second  (nonlinear)  order  in  ej,  fj: 

n  n  n 

Ci  =  X  X  Aijk  Q  fk  +  X  X  Bijkejek  +  X  X  Eijktjek 

j  k  j;^i„  j  k  j  k 

n  n  n  n  n  n 

fj  =  X  X  Cijk  Q  fk  +  X  X  Dijkejek  +  X  X  FijkQek 

j  k  j  k  j  k 

Here,  the  quantities  ei,  fi  on  the  lh.s.  are  regarded  as  dependent  variables,  while  ej,  fr  on  the  r.h,s.  are 
taken  to  be  independent  variables.  The  rule  for  the  selection  of  dependent  variables  is  this: 

(a)  if  an  independent  power  variable  is  fj,  the  selected  dependent  variable  is  ej 

(b)  if  an  independent  power  variable  is  ej,  the  selected  dependent  variable  is  fj. 

The  ei,  fi  equations  noted  above  are  the  canonical  equations  of  classic  theory.  In  drawing  multiport 
transducers,  the  symbol  TD  is  taken  to  be  the  bond  graph  structure  of  interlinked  junctions,  i.e.  the  interior 
system. 


Bulletin  #4  -  The  Electrostatic  Transducer  (A  2-Port  Model) 

(Ref.:  F.  Hunt,  “Electroacoustics”) 

The  simplest  embodiment  of  the  electrostatic  transducer  is  a  parallel  plate  capacitor  having  a  fixed 
separation  d,  a  fixed  charge  qo,  a  fixed  voltage  Eo,  and  an  area  S.  When  a  squeezing  force  is  applied  fluctuat¬ 
ing  in  magnitude  and  time,  a  flucutating  charge  q  is  induced  by  the  changing  separation  distance  x,  thus 
generating  a  fluctuating  voltage  Ee.  The  interaction  between  mechanical  and  electrical  variables  is  found  to  be 
nonlinear. 

The  Elementary  Bondgraph 

The  electrostatic  transducer  under  time  varying  input  is  modeled  by  transduction  T  with  power  variables 
F^,  X,  and  E^,  q: 

i  T  i;  F,  =  F,(q);  E,  =  E,(qx) 

A  q 
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Symbol  F,  is  the  force  (units:  Newtons)  between  the  plates  resulting  from  fluctuating  plate  charge  q. 
Symbol  is  the  electrical  potential  (units:  Volts)  resulting  from  the  combined  effects  of  fluctuating  charge  q 
and  fluctuating  distance  x  between  the  plates.  The  nonlinear  functional  relations  between  and  q,  and 
between  F^  and  q  x  are: 


P  _  1  [Jqdtf.  p  _  IJqxdtdt 

2  e„S  ’  ^  eoS 

Symbol  is  the  permittivity  of  free  space  (units:  coulombs/voltmeter).  From  these  choices  of  units, 
transduction  f  is  derived  to  have  the  units:  Volt  meter/Newton  Coulomb. 

A  simple  acausal  bond  graph  of  a  generalized  electrostatic  transducer  may  now  be  sketched: 


^ ^  source  Eg 


This  graph  may  be  read  as  follows:  an  applied  exterior  source  force  F^  overcomes  the  reactive  mechanical 
impedance  of  the  transducer  (namely  I^,  C^,  R  J  leaving  a  residual  force  F^  due  to  the  attraction  (or  repulsion) 
between  the  plates  resulting  from  the  effects  of  fluctuating  charge.  The  direction  of  the  bond  of  F^  is  left  for 
further  discussion.  The  transduced  voltage  E^,  q  to  satisfy  a  reactive  load  Zjj.  The  transduction  symbol  T  is 
bonded  to  two  1-junctions,  hence  it  simulates  a  0-junction,  i.e.,  it  acts  as  a  “coupler. 

Biased  Operation 

In  practice  electrostatic  transducers  operate  with  permanent  electrostatic  bias,  furnished  in  the  simplest 
structures,  by  an  electric  battery.  To  illustrate  biasing,  we  choose  the  case  of  an  elementary  capacitor  mode  of 
a  flexible  diaphragm,  area  S,  electroded  on  its  inner  surface,  set  at  a  distance  d  from  a  perforated  metal  screen 
which  is  at  ground  potential.  When  a  d.c.  bias  voltage  is  applied  to  the  capacitor  in  the  form  of  E^  (volts),  a 
charge  q^,  will  appear  on  the  plates,  and  simultaneously  a  steady  diaphragm  displacement  xo  from  equi^rium 
will  occur.  For  a  given  mechanical  compliance  Cm  of  the  diaphragm  there  will  be  an  effective  mechamcal 
force  XyC^  tending  to  restore  the  diaphragm  to  its  initial  shape.  In  the  absence  of  applied  (=  extenial)  me¬ 
chanical  forces,  this  restoration  force  will  be  balanced  by  the  electrically  generated  force  of  attraction. 


_  x«  Qo^  .  \r  _  q<i  (d+Xp) .  ^  ^ 

^0-  cc  ’  '^0 


CqS 

C„  ■  2£«S  ’  ■  ”  EoS  ’  d+x, 

These  equations  lead  to  an  implicit  equation  for  x^  of  third  degree. 


.  X  =  c  V^ 

■^0  ''o 


EqS 

2(d+x„)" 


Values  of  x^,  can  be  obtained  by  graphical  procedures. 

Two  modes  of  biased  operation  are  used  in  practice 

(a)  transmitter  mode,  in  which  an  electric  source  generates,  through  transduction,  an  acoustic  field. 
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(b)  receiver  mode,  in  which  an  acoustic  field,  acting  as  a  mechanical  source  generates  an  electrical  signal. 

Restricting  attention  to  the  receiver  mode,  it  is  assumed  that  a  fluctuating  acoustic  pressure  acting  over 
area  S  of  the  diaphragm  is  equivalent  to  a  concentrated  fluctuating  force  F,.  Such  a  force  generates  a  time 
fluctuating  displacement  X,,  together  with  a  time  fluctuating  charge  q,.  At  any  time  then,  total  displacement  X 
and  total  charge  q  are 


X  =  d  +  Xg  +  X,;  q  =  +  q, 

Expanding  q-  and  qx,  and  neglecting  second  order  terms,  one  finds  the  first  order  terms  in  fluctuating  q^ 
X|  to  be  2  q^qj  and  q^^x,.  The  forces  and  voltages  of  electromechanical  coupling,  to  first  order,  namely 
[read,  mechanical  force  due  to  a  fluctuating  ql]  and  [read,  electrical  voltage  due  to  fluctuating  xj  are 

p  _ qoQi  _  “1|  ,w  _  qoXj  _  Vj  ,  p,  _  €oS 

^"*^W"eoS  EoS  qo 

The  units  of  are  nW  or  C/N  (meters/volt  or  coulomb/Newton).  These  expressions  allow  one  to  form 
acausal  equations  by  use  of  the  bond  graph.  At  the  1 -junction  of  the  mechanical  side  of  summation  of  forces 
on  the  diaphragm  leads  to  F,  =  v,  ±  The  ambiguity  in  the  sign  of  Fme(l)  is  resolved  by  noting  that 
the  time  fluctuating  charge  ql  acts  as  a  source  so  that  the  force  balance  is  F,  +  F^^j^  =  Vj,  from  which 

Ft  =  ^  +Z^Xi  =  ^  +  Z^vi  =  T„e(I)  +  ZmVi 

Here,  T_^  (which  is  l/j(0  is  the  transduction  coefficient.  This  choice  of  sign  is  justified  on  physical 
grounds.  When  Fj  increases,  closing  the  gap,  the  electrostatic  force  of  attraction  also  increases;  similarly, 
when  F,  decreases,  opening  the  gap  (between  diaphragm  and  screen),  the  electrostatic  force  decreases.  Thus, 
the  fluctuation  of  charge  appears  on  the  mechanical  side  as  a  force  source  adding  to  the  applied  force  Fj. 

On  the  electrical  side  of  transduction  the  electrical  effect  of  gap  motion  in  the  form  of  displacement  Xj  is 
^emd)  “  =  X,/C^^.  This  is  the  potential  generated  by  the  fluctuation  X,.  It’s  direction  is  awav  from  the 

1-junction  labelled  q,  that  is,  it  constitutes  a  negative  source.  In  absence  of  an  applied  electric  potential  at  the 
electrical  terminals,  summation  at  the  1-junction  q  leads  to. 


Z  qi  =  Vnjg(i),  Z  =  Zg  +  Z^,  Zj(  =  electrical  load 


In  terms  of  first  order  fluctuating  variables,  this  summation  of  potentials  is, 

0  =  Zqi+-^ 

Assume  now  that  fluctuation  in  time  is  sinusoidal.  By  elimination  of  q,  one  finds 


Fi  =  X,[Z„ 


(i® 


] 


In  absence  of  an  electrical  load,  at  low  frequencies, 

Z  =  1  •  7^  —  1  • 
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Thus  as  CO  0, 


Fi 

00 ->0 


^nPo  1  _ 

^cm 


Here  the  quantity  c  C  /C  -  describes  the  coupling  between  mechanical  and  electrical  ‘compliances’  at 

~  •'me  cm  i  /  j\ 

very  low  frequency.  It  is  designated  as  k^  the  coefficient  of  electromechanical  coupling  (squared), 


.  v2  —  •  C'  —  ^0 

Q2„  ’"-o-  dtxo’'’'®"'"  q„  Eo 

These  formulas  show  to  be  a  function  of  E^,  q^,,  X^,  and  of  diaphragm  compliance  C^.  From  energy 
considerations  there  is  a  limit  on  the  magnitude  of  k^ 

>  1 

This  is  also  a  requirement  for  stability  of  the  transducer.  For  given  the  achievable  bias  displacement 
is  limited, 


1  Eq^S 

C„  "  (dtxo)’ 

It  can  also  be  deduced  from  these  formulas  that  the  low  frequency  electrical  capacitance  C',  also  called 
the  ‘free-running  capacitance’  is 

Cj  =  ,  where  C,,  is  the  ‘blocked  capacitance.’ 

Equivalent  Circuit:  Ideal  Transformer:  Turns  Ratio 

Although  every  equivalent  circuit  may  be  transcribed  into  a  bond  graph,  a  bond  graph  is  not,  per  se,  an 
equivalent  circuit.  Steps  required  to  form  a  2-port  equivalent  circuit  of  an  electrostatic  transducer  are  outlined 
next 

(1)  Let  Ej,  F,  be  the  terminal  electric  potential  and  mechanical  force  respectively  of  the  2-port.  On  the 
electrical  side  the  blocked  electrical  capacitance  €„  is  moved  to  ‘shunt’  position,  that  is,  from  a  1-junction  to  a 
0-junction: 


The  purpose  of  this  move  is  to  make  C^,  a  ‘coupler’  between  I,  and  Vj. 

(2)  A  ratio  ())  is  formed  of  the  transduction  coefficient  T^^  and  Zg=  l/j(0  C^,, 


<l>  = 


1/jCO  Oem  Op 

l^Cd  Co  Ccm 


c 

(units:  —  or 
m 


li 

V 


) 
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This  quantity  is  the  “terms  ratio”  of  an  ideal  electromechanical  transformer  which  relates  the  independent 
variables  I,,  v,  of  the  canonical  equations  to  the  dependent  variables  F,,  E,: 

E,^  =  F,:  I,  =  (t»v,;  E,/I,=F,/v,(t.Me.  = 

In  terms  of  the  bias  quantities  E^,,  X^,  q^, 

.  _  EqCq  _  go 

^  di  x„  dt  x„ 

By  shifting  to  the  shunt  position,  the  mechanical  compliance  changes  from  cm  to  C^,  a  result  of 
making  the  blocked  capacitance  the  coupling  element. 

The  2-port  equivalent  circuit  is  now  complete:  its  bond  graph  is, 


The  arrow  directions  are  conventional  for  a  2-port  electromechanical  quivalent  circuit. 


Bulletin  #5  -  Mathematical  Modeling  of  Electromechanical  Transducers 
Based  on  Electrical  Four-Pole  Theory 

(Ref.:  “Grundlagen  DerTechnischen  Akustik,”  W.  Reichhardt,  Leipzig,  1968) 

In  electrical  4-pole  theory  the  power  variables  are  Ej,  ij,  E^,  ij.  Chosen  as  independent  variables  are  E^,  i^, 
which  implies  the  functional  relations, 

E,  =  <P,(E2,i2) 
ii  =  <p2<E2-i2) 

To  a  first  linear  approximation,  one  finds, 


Ej  —  AJ2E2 

'e.' 

^2 

’a„ 

^12 

“^21  ^2 '^'^22  ^2 

_ii  _ 

=  A 

_Ea_ 

;  A  = 

A21 

^22 

In  general  all  power  variables  are  complex  phasors  in  steady  state  operation,  requiring  (then)  all  coeffi¬ 
cients  Aij  to  be  complex  numbers. 

There  are  limits  to  the  components  Aij.  These  are  set  by  the  assumption  that  the  transfer  of  power 
through  the  4-pole  is  lossless.  It  is  then  seen  that, 

|(E,  i>E;i,)-i(]^2i2+E^2)  =  0 
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From  this  requirement  one  deduces  that, 


(a)  (A„A,>A*A,,)i,i:  =0 

(b)  (A,3A,*+A*A^)E,E:=0 

(c)  (A„A2*+A*A,,-l)i,E*  =0 

(d)  (A„A>A*  A^-l)i*E,  =0 

Two  types  of  transducers  can  be  modeled  by  the  El,  il,  equations  noted  above: 

(1)  El  oc  i2  for  the  (short  circuit  condition)  =  0.  In  from  (a)  above  A^,,  is  purely  imaginary  (written 

j  /^).  A„  (at  our  option)  is  made  purely  real  which  makes  A,^  purely  imaginaiy  (written  A^),  as  seen  from 
(b).  These  assignments  are  inserted  into  (c)  resulting  in 

A/iAi-jAfJjAii  =  |A|  =  +  l, 

(2)  In  the  second  type  E, «  E^  for  the  (open  circuited)  conduct  i^  s  0.  In  this  case  Aj^  is  purely  real 
(written  A^, )  because  E,  is  assumed  purely  real.  From  (b)  above  it  is  seen  that  A^  is  purely  imaginary  (written 
A^).  Aj,  (4  our  option)  is  made  purely  real,  which  from  (a)  above  makes  Aj,  purely  imaginary  (written  j 
A^).  These  assignments  inserted  in  (d)  result  in  the  formula 

Ai4A,VjA/|jA^-l=0 

or 

jAn  jA^-  Aj4  A^,  =  |A|  =  -1 

This  expression  for  the  determinant  of  A,  as  well  as  the  one  in  (1)  above  can  serve  to  characterize  electro¬ 
mechanical  transducers. 


Four-Pole  Modeling  of  Electromechanical  Transduction 

The  transduction  component  of  a  transducer  can  be  modeled  as  a  4-pole.  In  bond  graph  rotation,  the 
power  variables  Em,  v,  Fe  appear  as  bonds. 


.Melt  — 

i  ^  Fe 

E^  =  electric  field  coupled  to  the  mechanical  circuit 
F^  =  mechanical  force  coupled  to  the  electrical  circuit 
i  =  current 

V  =  mechanical  velocity 

Relations  between  variables  are  based  on  v,  Fe,  chosen  as  independent  variables,  and  Em,  i,  as  dependent 
variables. 


E  =A,,v+A,,F, 

m  II  12  e 


i  —  Ajj  V  +  A22  F^ 


Two  types  of  transduction  may  be  deduced. 
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( 1 )  In  the  transfer  term  Aj,  v  one  sets  v  =  0.  Then  Em  =  Aj ,  Fe.  The  determinant  | A|  =  - 1 .  This  type  of 
transduction  characterizes  the  magnetic  field  transducer 

(2)  In  the  transfer  term  A,,  v  one  sets  v  s  0.  Then  Em  =  Aj,  Fe.  The  determinant  | A|  =  - 1 .  This  type  of 
transduction  characterizes  the  electric  field  transducer. 

The  transduction  component  is  one  link  in  a  chain  (=  cascade)  of  transducer  components,  every  link  of 
which  is  itself  a  4-pole. 

Four-Pole  Low  Frequency  Canonical  Equations  and  Corresponding  Bond  Graphs 
Magnetic  Field  Transducers 

A  low  frequency  bond  graph  has  the  general  structure  given  by 

—  1  — T  —  0  —  (=  “one-tee-zero”) 


In  detail. 


E:Se 


U 

I 

T.  1  -f=!^T 


Se;F 


(Arrows  show  transmission  from  electrical  field  to  mechanical  field.) 

Explanation 

(1)  The  source  Se  at  the  left  is  a  voltage  (=  effort)  source  and  must  therefore  feed  into  a  1 -junction. 

(2)  The  source  Sf  at  the  right  is  a  force  (=  flow)  source  and  must  therefore  be  bonded  to  a  0-junction. 
This  is  the  result  of  representing  the  mechanical  side  by  its  dual  form. 

(3)  The  blocked  coil  inductance  Lb  is  bonded  to  the  1-junction  (an  effort  divider)  to  represent  the 
physical  requirement  that  for  the  condition  (o^— >0  the  voltage  Ej^->0,  leaving  source  E  to  be  equal 
to  Em.  If  Lf  was  attached  to  a  0-junction  the  source  would  be  short-circuited  and  no  d.c.  current 
could  flow. 

(4)  The  mechanical  compliance  Q  is  bonded  to  the  common  velocity  (dual  form)  0-junction  to 
represent  the  physical  requirement  that  for  the  condition  ©^0  the  force  F^  becomes  infinite 

(=  “open  circuit”)  and  the  velocity  (i.e.  the  a.c.  component)  becomes  zero  (=  clamped  condition). 
Bonding  Q  to  a  1-junction  (in  dual  form)  would  mean  v^->0  as  ©^0  so  that  v^-»n  leading  to  a 
finite  velocity  at  zero  frequency,  which  is  unphysical. 

(5)  Writing  for  brevity  the  transduction  coefficient  as  X,  one  may  codify  the  transduction  equations 
as. 


i  =  XF, 
F  = 

m  X 
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The  bond  graph  equations  are  found  by  summing  voltage  at  the  1 -junction  and  force  at  the  0-juntion, 


E  =  jwLGi  +  -^'' 

i  =  X  +  F] 

L  jw  Cn 

Substitution  of  i  into  the  E  equation  leads  to  the  canonical  set  in  v,  F,  whose  matrix  is, 

A  =  /A,-L  jAii\ 

The  determinant  |A|  =  +1,  as  required. 

Electric  Field  Transducers 

A  low  frequency  bond  graph  has  the  general  structure  given  by, 

—  0  —  T  —  0  —  (=  “zero-tee-zero”) 


In  detail. 


Cb 

c 

ME 


Q 

c 

V>|F, 


■0 


E  T  ''  n  '' 
■T  n'V 


im 


Fe 


It  is  seen  that  the  electrical  field  is  in  “direct  form”  and  the  mechanical  field  is  in  “dual  form.”  The  choice  of 
dual  form  comes  from  the  choice  of  v,  F  as  independent  variables. 

Explanation 

(1)  The  blocked  capacity  is  bonded  to  a  0-junction  to  represent  the  physical  requirement  that 
for  the  condition  m-^O  the  current  i^  vanishes  and  the  “supply  current”  i  becomes  the  current 
**  coupled  to  the  mechanical  circuit. 

(2)  The  source  driving  the  blocked  capacity  must  be  a  ‘current  source’ S,.  Ordinarily,  this  type  of 
source  is  furnished  by  a  ‘voltage  source’ bonded  to  a  1-junction  which  is  loaded  with  a 
source  impedance  Z^: 

Zs 

ilE. 

E«  E 

Sf  =  Se  1  ■  ^  (=  equivalent  current  source) 
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(3)  Writing  for  brevity  the  transduction  coefficient  as  Y,  one  may  codify  the  transduction  equa¬ 
tions  as, 


E=YF,  V  1  m 

.  1  > 


The  bond  graph  equations  are  found  by  summing  currents  at  the  0-junction  and  forces  at  the  (dual) 
0-juntion: 


i=jw(;E+^'' 

E  =  YFe  =  Yr:-^  +  F] 

•-JW  Cm  J 

Substitution  of  E  into  the  i  equation  leads  to  the  canonical  set  in  v,  F,  whose  matrix  is, 


A  = 


iAii  =  _ X _ •  AX  =Y 

"  jWCn’  “ 


The  determinant  |A|  =  -1,  as  required. 

Conclusions 

Modeling  electromechanical  transducers  on  the  basis  of  electrical  4-pole  theory  involving  the  power 
variables  E,  i,  v,  F,  leads  to  a  set  of  canonical  equations  in  which  choices  of  independent  variables  must  be 
made.  For  various  reasons  several  authors  choose  v  and  F,  with  the  additional  choice  of  ‘v  across’  and  ‘F 
through.’  The  mechanical  branch  is  then  represented  in  dual  form.  Other  authors  take  i,v  to  be  independent 
variables,  both  in  direct  form,  namely  ‘i  through’,  ‘v  through.’  The  resulting  canonical  set,  while  different 
from  the  set  based  on  4-pole  theory,  can  be  converted  to  it  by  interchange  of  variables. 

Basic  transducer  types  are  simply  summarized: 

Magnetic  field;  —  1  —  T  —  0 — ;  driven  by  a  “voltage  source” 

Electric  field;  —  0  —  T  —  0  — ;  driven  by  a  “current  source” 

Coefficient  of  Electromechanical  Coupling  and  Its  Effect  on  Low  Frequency  Properties  of  Inductance. 
Capacitance  and  Compliance.Case  I.  F.  V  at  right  angles  to  E.  D  or  B.  H. 

The  bond  graph  that  describes  the  low  frequency  behavior  of  a  magnetic  field  transducer  calls  for  a 
blocked  inductance  Lb  bonded  to  a  1-junction  and  a  mechanical  compliance  bonded  to  a  0-junction  (dual 
form).  Because  of  electromechanical  coupling  the  magnitude  of  the  compliance  depends  on  the  nature  of  the 
electrical  drive.  Since  magnetic  field  transducers  are  driven  by  a  voltage  source  a  purely  (uncoupled)  me¬ 
chanical  magnitude  of  compliance  can  be  measured  by  opening  the  electric  input  circuit.  The  open  circuit 
compliance  is  then  written  a  Cj(.  When  the  circuit  is  closed  and  a  current  i  flows  through  the  (coil)  inductance 
Ljj  the  open  circuit  compliance  transferred  to  the  electrical  branch  appears  as  coupled  inductance 
bonded  to  the  1-junction.  The  sum  where  is  the  “free-running”  inductance.  By  defining  the 

coupling  coefficient  =  Ly(L^  +  L^^),  it  is  readily  deduced  that  the  measured  free-running  inductance 
Lf  =  Lyi-k^  Further  more  it  may  readily  be  demonstrated  that  the  measured  compliance  C^.  on  short  circuit 
drive  is  related  to  the  open  circuited  compliance  Q  by  the  formula,  C^=  Q  (1-k^).  It  also  may  readily  be  seen 
that  because  k^  <  1  the  short  circuit  compliance  is  less  than  the  open  circuit  compliance. 
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Similar  considerations  can  be  applied  to  electric  field  transducers.  Since  they  are  driven  by  a  current 
source  the  mechanical  compliance  q  is  measured  at  electrical  short  circuit  thereby  bypassing  the  blocked 
capacity  C ..  When  a  drive  voltage  is  restored  the  mechanical  compliance  decreases  to  the  free-running 
magnitude  Q  =  q  (l-k^).  Correspondingly,  the  “free-running”  electrical  capacity  C,  is  increased  to 
q=Cb/l-k\ 

Case  II.  F.  v  parallel  to  E.  D  or  B.  H 

In  this  case  negative  compliances  appear  (see  Bulletin  on  the  electrostatic  transducer).  Summarizing  the 
changes  to  be  made  to  Case  I: 

Electric  field  transducers 

short  circuit  compliance  C^  is  replaced  by  open  circuit  compliance  C;  plus  a  negative  compliance  - 
Ci/k^;  i.e. 


-^0-^  is  replaced  by 


Case  I  Case  11 


Magnetic  Field  Transducers 

open  circuit  compliance  C/  is  replaced  by  short  circuit  compliance  q  plus  a  negative  compliance  - 
q/k^;  i.e. 


-|-0^  is  replaced  by 
Case  I  Case  II 


The  rationale  for  these  changes  is  explained  by, 
A.  Lenk,  Acoustica  (1958)  Nr.  3,  page  159 


Rnllfttin  #6  -  Moving  Armature  Transducer 
(Ref.:  “Electric  Circuits,”  MIT,  J.  Wiley  &  Sons,  Inc.,  New  York,  1943) 


The  Moving  Armature  Transducer  (shown  in  the 
simplified  sketch  below)  is  based  in  operation  on  energy 
transfer  between  electrical  and  mechanical  ports  associ¬ 
ated  with  time  variation  of  reluctance  in  a  magnetic 
circuit.  As  shown  it  has  three  ports:  (1)  a  field  winding 
of  electrical  wire,  q  turns,  around  an  iron  yoke.  Its 
power  variables  are  potential  difference  v,  (units:  volt) 
and  current  i,  (units:  amp.)  (2)  an  armature  winding  of 
electrical  wire,  N  turns,  around  the  pivoted  iron  armature. 


21 


Its  power  variables  are  v  (units:  volt)  and  i  (units:  amp).  (3)  a  mechanical  port  near  the  gap  of  the  magnetic 
circuit.  Its  power  variables  are  force  F  (units:  newton)  and  velocity  V  (units:  meter/sec). 

In  operation  the  gap  energy  Wm  for  a  magnetic  feed  strength  H  (units:  amp/meter)  and  magnetic  flux 
density  B  (units:  Vs/m-)  is  (HB/Stt)  Sx  =  W^,  where  Sx  is  the  gap  volume.  When  a  force  F  is  applied  at  the 
mechanical  port,  it  changes  the  gap  with  x.  The  gap  magnetic  energy  changes  in  the  amount  dW^dx  =  F. 
Allowing  H  =  B/)xo  in  the  gap  it  is  seen  that  F  =  SBV8p  po. 

At  the  electrical  port  of  the  field  winding  energy  can  only  be  extracted  in  the  amount  i^  d  (Nj(p)/dt.  The 
field  current  i^  is  nearly  constant  and  the  flux  (p  (units:  Vs),  vaiying  around  a  constant  mean  value,  is  also 
nearly  independent  of  time  to  the  approximation  used  in  this  analysis.  Thus  over  time  the  field  energy  W^.  at 
the  field  port  does  not  change. 

At  the  electrical  port  of  the  armature  winding  the  rate  of  energy  extraction  when  armature  current  i  is 
active  is  dWa/dt  =  Nidtp/dt.  Since  current  i  is  variable  due  to  variable  F,  and  flux  tp  in  the  gap  varies  with 
changing  gap  5x,  then  energy  may  be  transferred  at  the  armature  winding.  The  total  flux  in  the  gap  is  the  sum 
of  a  constant  component  due  to  i^  and  a  variable  component  Ni 


4a:(N^^-fN.)  Sp„  ^  4?!  (1 -f  Nf ) 

(X„-^AX) 


If  Ni/Ni^  and  8x/Xo  are  small  relative  to  their  squares,  when  the  division  is  made  and  second  order  terms 
rejected  the  flux  in  the  gap  is  approximated  by 

4,iSm<,NA  R 

X„  N,1,  X„ 

Since  this  flux  varies  with  time  and  displacement  there  is  a  voltage  v  induced  in  the  armature  winding, 

v  =  N—  =  L,  —  =  B,  =  —  V  (v  =  velocity) 
dt  dt  o 

where  SpoN  .  g  _  47t  j^Nfif  (units:  Vs/m^) 

o  Xo 

It  is  seen  that  is  the  self-inductance  of  the  armature  winding.  Thus  the  induced  armature  voltage  has 
two  components:  the  first  in  di/dt  is  the  voltage  of  self  induction,  and  the  second  an  e.m.f.  generated  by 
mechanical  motion. 

It  was  noted  earlier  that  the  applied  force  in  the  magnetic  circuit  air  gap  is  F  =  S  BV8jt  po,  where  B  is  the 
flux  density.  Noting  that  B  =  rp/S  and  rejecting  (Ni/N^i^)^,  (Ax/xo)^,  (Ni/N^^)  (Ax/xo),  as  terms  of  second  order, 
it  is  seen  that  to  first  order, 

P  _  2jt  Spo  (Nfif)^  27t  Ax  +  2Ni . . . ) 

O 
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Each  term  in  this  3-term  form  has  significance.  (1)  the  first  term  states  that  there  is  a  component  of  force 
due  to  current  i^  of  fixed  magnitude  tending  to  close  the  gap.  (2)  the  second  term  states  that  there  is  a  negative 
component  of  force  proportional  to  Ax  which  tends  to  increase  the  displacement,  that  is  the  application  of  this 
force  widens  the  gap.  (3)  the  third  term  states  that  there  is  a  component  of  force  that  depends  on  armature 
current  i,  that  is,  when  time  varying  current  is  is  impressed  in  the  armature  ciruit  it  causes  the  armature  to 
move  in  unison  with  it,  and  when  the  armature  is  moved  by  a  mechanical  force  an  e.m.f.  in  the  electrical  port 
is  generated  as  noted  above. 

Construction  of  the  Bond  Graph  (Ref.;  Reichhardt,  see  above). 

Since  the  field  current  i^  is  virtually  constant  and  the  armature  current  i  is  variable  with  time  some 
analysts  combine  their  respective  ports  into  one  electrical  port,  the  power  variables  of  which  are  u  (units:  volt) 
and  i  (units;  amp)  where  i  ^  =  i^  +  i.  We  adopt  this  procedure  here.  Now  the  magnetic  flux  due  to  ip^  is  bPl/ 
Rp^,  where  Rp^  is  the  reluctance  of  the  magnetic  circuit  (units:  C/VS^).  The  input  admittance  of  the  combined 
el^trical  port  is  Rj^/j  where  is  the  total  turns  of  the  field  wire  (units:  none).  The  first  element  of  the 
bond  graph  is  a  current  divider  (i.e.  a  0-junction) 


Lpe 

ipe 

^ - 0- 

Iges  u 


Ipc  — 


U  Rc 


jwN/ 


r 


The  second  element  of  the  graph  involves  the  armature  current  i.  At  the  gap  the  magnetic  flux  Ogap  = 
Nip.oS/xo  (units:  VS).  The  drop  in  voltage  across  the  gap  is  jco  NOgap  or  Joj  ^i,  ^  poS/Xo.  This  is  the 
input  impedance  of  the  current  i.  It’s  graph  is  a  1-junction: 

Lb 

(u-uj 

U  J  “w  •  T  • 

- 1 — ^  u-Uw=JwLbi 

i 

The  third  element  in  the  graph  is  the  transduction.  In  electromagnetic  circuits  it  is  a  gyrator,  defined  as  the 
‘gyration’  of  electrical  voltage  (u)  into  mechanical  velocity  v  (units:  m/s),  and  the  ‘gyration’  of  electrical 
current  into  mechanical  force  F  (units:  Newton).  Its  graph  is 


rw 


u^  =  X“"v; 


Fw  =  X*i 


Here,  U^,  is  an  ‘across’  variable  and  is  a  ‘flow’  variable. 
From  earlier  discussion. 


X*  =  Bf  ^  where 


„  _  4jiNfifPo 
Bf - ^ - 


V 

(units;  — ^) 
m^ 


and,  X  has  the  units  of  Vs/m. 
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The  fourth  element  in  the  graph  is  the  mechanical  circuit.  In  this  circuit  is  armature  described  dynamically 
by  lumped  parameters  of  mass  M  stiffness  K  mechanical  resistance  and  radiation  resistance  R^j,.  These 
are  all  attached  to  a  flow  (i.e.  force)  divider,  which  is  a  0-junction: 


1 


1 

Rrad 


Here,  F^,  is  the  force  of  the  negative  stiffness  discussed  earlier. 
The  Completed  Graph 
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Application  to  the  Telephone 

In  the  telephone,  the  human  voice  supplies  a  variable  pressure  to  a  metal  disphragm.  This  is  the  armature. 
It  vibrates  with  a  velocity  which  is  transduced  into  an  electrical  current,  which  is  transmitted  to  a  distant 
receiver.  At  the  receiver  a  voltage  is  generated  upon  reception  of  current.  This  voltage  is  transduced  into  a 
mechanical  velocity  of  a  diaphragm  which  radiates  sound,  completely  reproducing  the  message  sent  by  the 
transmitter. 
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